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Abstract. Spin wave excitations and stability of the (0, tt) ordered SDW state are 
investigated within the minimal two-band model for iron pnictides including a Hund's 
coupling term. The SDW state is shown to be stable in two distinct doping regimes 
— finite hole doping in the lower SDW band for small second neighbour hoppings, 
and small electron doping in the upper SDW band for comparable first and second 
neighbour hoppings. In both cases, Hund's coupling strongly stabilizes the SDW state 
due to generation of additional ferromagnetic spin couplings involving the inter-orbital 
part of the particle- hole propagator. Spin wave energies for the two-band model are 
very similar to the one-band t-t' Hubbard model results obtained earlier, and are in 
agreement with inelastic neutron scattering studies of iron pnictides. 
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1. Introduction 



Single- crystal neutron scattering studies of iron pnictides have indicated a commensurate 
magnetic ordering of iron moments ordered ferromagnetically in the b direction and 
antiferromagnetically in the a and c directions [l]. Inelastic neutron scattering 
measurements in AFe2As2 (A = Ca, Ba, Sr) yield well-defined spin-wave excitations up 
to the zone boundary on an energy scale ~ 200meV [2l[3lll]. The realization of a (0, tt, vr) 
ordered SDW state has opened the possibility of observing phenomena in this class 
of compounds which are characteristically associated with both the antiferromagnetic 
(AF) state such as quantum spin fluctuations, hole/electron motion in AF background, 
spin-fluctuation mediated pairing as well as the metallic ferromagnetic (F) state such 
as carrier-induced spin interactions, correlation-induced spin-charge coupling, and non- 
quasiparticle states. 

A minimal two-band model consisting of two degenerate orbitals dxz and dy^ 
per Fe ion on a two dimensional square lattice [5] has been widely studied recently 
O m El El [ini [n] in order to understand the magnetic ordering and excitations in doped 
iron pnictides. Ab initio calculations [12] using local density approximation suggest that 
the Fermi surface is determined by bands having mostly dxz and dyz character, as also 
indicated by polarized angle-resolved photoemission spectroscopy (ARPES) [13]. The 
hybridization of Fe 3d orbitals with themselves as well as through the As 3p orbitals 
lying above and below the square plaquettes formed by the Fe atoms leads to effective 
hopping parameters of the two-orbital model as shown schematically in Fig. 2. 

However, the role of inter-orbital exchange interaction (Hund's coupling) on 
spin wave energies and the stability of the (0, n) ordered SDW state have not been 
investigated within the minimal two-band model. This study would also allow for 
comparison with the one-band model results obtained earlier within the t-t' Hubbard 
model [T^ [T5]. The F-AF state in two and three dimensions was shown to be stabilized 
by the AF spin couplings generated by the second neighbour hopping t' as well as 
the carrier-induced F spin couplings as in metallic ferromagnets, which are strongly 
enhanced by the t' induced asymmetric and peaked electronic spectral distribution due 
to band dispersion saddle points. 



2. (0, tt) ordered SDW state of the two-band model 

Hund's coupling is known to stabilize metallic ferromagnetism in orbitally degenerate 
systems [16]. Recent investigations of the correlated motion of electrons in multiband 
ferromagnets [17| have yielded an effective quantum expansion parameter ^p^l^^~^p^2 
in terms of the orbital degeneracy A/", intrasite Coulomb interaction U, and Hund's 
coupling J, which provides an effective measure of the correlation induced quantum 
corrections which reduce the spin stiffness and spin wave energies and thus destabilize 
the ferromagnetic state. The strong suppression of this quantum parameter by Hund's 
coupling, especially for large A/", provides insight into the stabilization of metallic 
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Figure 1. The effective hopping parameters of the two-band model involving the two 
orbitals dxz and dyz per Fe ion on a square lattice (from Ref. [5]). The different cr and 
TT orbital overlaps result in anisotropic intra-orbital hoppings ti and t2 with opposite 
sign in different directions. Also included are second neighbour intra-orbital (t^) and 
inter-orbital {t^) hoppings. 



ferromagnetism in orbitally degenerate systems. As these quantum corrections involving 
self energy and vertex corrections will affect the ferromagnetic spin couplings, the role of 
Hund's coupling on magnetic excitations in doped iron pnictides is therefore of interest. 

The inter-orbital Hund's coupling term can be conveniently included on an equal 
footing with the intra-orbital Hubbard interaction term in a general multi-orbital 
correlated electron model: 



where the interaction matrix elements f/^j/ = ioi jj, = p and U^^p = 2J for n ^ v 
refer to the intra-orbital and inter-orbital Coulomb interaction terms, respectively. In 
the following, we will specifically consider a two-orbital model on a square lattice with 
reference to the dxz and dy^ orbitals of interest for iron pnictides, with the four hopping 
terms ti - t^ as shown in Fig. 1. The inter-orbital density interaction term Vrii^Uiy has 
been dropped as it does not play any role in the magnetism up to the random phase 
approximation (RPA) considered here. 

We consider the (0,7r) ordered SDW state of the above two-band model, with F 
and AF spin orderings along the x and y directions, respectively. For the single-band 
model, this SDW state can be conveniently represented in a two-sublattice basis [I^ITSj . 
Extending this approach to a composite two-orbital (a two-sublattice (A B) basis, 
the Hartree-Fock (HF) level Hamiltonian matrix in this composite basis (Aa A/3 Ba 



H 




(1) 
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B/3) assumes the form: 
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where a and refer to the dxz and dy^ orbitals, the band energies 



Ax 



.ly 



-2ti cos 
-2ti cos ky 



2x 



— 2t2 COS k^ 

— 2t2 cos /Cy 



-4^3 cos /cj; COS A;,, 



-4^4 sin A;^. sin fc,, 



corresponding to different hopping terms in different directions, and the self-consistently 
determined exchange fields: 

2Ap = Upmp + Jm„ (4) 

in terms of the sublattice magnetizations mQ and for the two orbitals. Quantum 
corrections to sublattice magnetization resulting from inter-band spectral weight transfer 
due to electron-magnon interaction should contribute to the substantially reduced 
magnetic moment observed in iron pnictides. 

It is instructive to consider two limiting cases which connect to the SDW state of 
the one-band model. When ti = ^2 = t, t^ = t', and ^4 = 0, the two orbitals in Eq. (2) 
get decoupled and are also exactly degenerate as ti = ^2, and Eq. (2) identically reduces 
to the SDW state Hamiltonian for the one-band t-t' Hubbard model 1151: 



-a 



A 



a 



A 



(5) 



where e^^'^ = —2t cos kx(y) and = — 4t' cos /c^: cos A;^,. 

In the more relevant limiting case ti = — ^2 = —t, t^ = —t', and t^ = 0, although the 
two orbitals in Eq. (2) again get decoupled, they are no longer degenerate as ti ^t^- For 
the /3 orbital, although the off-diagonal terms have exactly opposite sign as compared 
to Eq. (5) for the one-band model due to the ti, ^3 sign reversal, the energies and 
amplitudes remain unchanged due to energy band folding in the SDW state. For the a. 
orbital, the transformation kr^ ^ k^ ^ cancels the effect of the ti, t3 sign reversal on 
both diagonal and off-diagonal terms, and again renders the SDW state identical to Eq. 
(5) for the one-band model. 

The above equivalence also implies identical spin-wave dispersion as for the one- 
band model. In the following, we will investigate the effects of finite band mixing term 
^4 on the spin wave energies in the (0, vr) ordered SDW state. 
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Figure 2. The spin-wave dispersion obtained for the two-band model with Hund's 
couphng (J — C//4), without Hund's coupling (J = 0), and for the single-band case 
(^4 = 0). Here t3=t4=— 0.3, and the hopping energy scale \ti \ = 200nieV. 



3. Transverse spin fluctuations in the broken-symmetry state 

In the random phase approximation, the transverse spin fluctuation propagator for the 
two-band model includes both U and J ladders, and hence retains its usual form: 

[XRPA(q,-)]-i_[^]^^0(q,,)] (6) 

where the interaction matrix \U] includes the Hund's coupling term J as off-diagonal 
matrix elements, as given below (1). Spin wave energies were obtained approximately 
from poles of Eq. (6) in terms of the largest eigenvalue of real part of \U] The bare 
particle-hole propagator: 



k,Lm 



(7) 



is evaluated in the orbital-sublattice basis by integrating out the fermions in the (0,7r) 
ordered SDW state. Here E^^ and ^ko- are the eigenvalues and eigenvectors of the 
Hamiltonian matrix (2), the orbital-sublattice basis indices a,b run through 1 to 4, 
and l,m indicate the four eigenvalue branches. The superscripts +(— ) refer to particle 
(hole) energies above (below) the Fermi energy, and both inter-band and intra-band 
particle-hole terms are included. 

The bare propagator [x^] has a flnite imaginary part representing (low-energy) 
intra-band and (high-energy) inter-band particle-hole excitations, resulting in flnite spin 
wave damping and linewidth even at the RPA level, as in metallic antiferromagnets 
[T8] . The spin-charge coupling mechanism relevant for metallic ferromagnets such as 
manganites will also contribute due to decay of spin waves into longer wavelength modes 
accompanied with internal charge excitations [19] . Indeed, this has been suggested from 
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Figure 3. (a) Hole doping dependence of spin-wave dispersion for the two-band model 
with Hund's coupling included. Here t3=t4=— 0.7 and A = 2. (b) Without Hund's 
coupling (J = 0), the spin wave energies become negative over a larger part of the 
Brillouin zone. 



the absence of any steep increase in damping at higher energy indicative of a Stoner 
continuum [3]. From the observed high energy behaviour of spin wave damping ascribed 
to particle-hole excitations [1] , it has been inferred that the full excitation spectrum can 
not be understood in terms of the local moment picture. Weakly damped spin waves 
near the ordering wavevector have been obtained within multiband models from the 
imaginary part of the spin fluctuation propagator [Sni EH |22] . 

Fig. 2 shows the spin- wave dispersion in the (0,7r) state of the two-band model, 
showing the stabilization of the F-AF state for small Hund's coupling. For J = 0, the 
spin wave energy becomes negative for small near the M and F points, indicating 
instability with respect to long wavelength spin twisting modes in the ferromagnetic 
ordering direction. Here, the hopping terms are ti=— 1.0, ^2=1-0, and ^3=^4=— 0.3, 
the exchange field A=3.0, and the hole doping concentration x ~ 35%, for which the 
self-consistency condition (4) yields U -f J ~ 11. Also shown is the single-band result 
obtained by simply setting t^ = 0. We have set = 1 as the unit of energy scale, and 
the spin wave energies are shown for |ti|=200meV. 

Evidently, the electronic spectral function modification due to the band mixing term 
^4 reduces the F spin couplings and destabilizes the F-AF state. The additional spin 
couplings J'^[x^]'^f generated by Hund's coupling due to the inter-orbital component of 
the particle-hole propagator restores the stability of the F-AF state. 

Fig. 2 also shows that for small t^, the overall structure of the spin wave dispersion 
for the two-band model, with Hund's coupling included, is very close to that for the 
one-band model obtained earlier. [T5] with respect to anisotropic spin wave velocities, 
spin-wave dispersion, and energy scale. Significantly, the spin-wave dispersion clearly 
shows a maximum at (tt, tt) in agreement with neutron scattering experiments [21 |3l H]. 
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Figure 4. (a) Spin wave energies for the two-band model with same hopping 
parameters as in Ref. [5], showing strong enhancement (suppression) in the F (AF) 
ordering direction with increasing electron doping y. (b) Spin wave energies are 
strongly enhanced when a small Hund's coupling term is included. 

For larger second-neighbour hoppings t3=t4=— 0.7 and somewhat smaller exchange 
field A=2, the hole doping dependence of spin wave dispersion is shown in Fig. 3(a). 
The self-consistency condition yields U + J ^ 8 for x = 20%. While the peak spin wave 
energy in the AF direction is indeed enhanced, as expected from the stronger second- 
neighbour AF spin couplings generated, the F-AF state is evidently not robust with 
respect to fiuctuations in the F ordering direction. For x = 10%, the spin wave energy 
becomes negative near X and M', indicating instability with respect to zone boundary 
modes, whereas at higher hole doping x = 30%, the spin wave energy becomes negative 
near F and M, indicating instability with respect to long wavelength modes. 

For the relatively stable case at x = 20%, the spin wave energies with and without 
Hund's coupling are compared in Fig. 3(b). Again, Hund's coupling is seen to play a 
crucial role in stabilizing the F-AF state, as indicated by the negative spin wave energies 
obtained for J=0 over a larger part of the Brillouin zone. 

Finally, we consider the set of hopping parameters ti=— 1.0, t2=1.3, ^3=^4=— 0.85, 
which yields circular electron and hole pockets near the bottom of the upper band. [5] 
At half filling, the spin wave dispersion (Fig. 4(a), ?/ = 0) is very similar to that 
derived for the single band t-t' model [H], for which the spin wave energy vanishes at 
the zone boundary in the ferromagnetic direction due to absence of any ferromagnetic 
spin coupling in the insulating state. 

For small electron doping in the upper SDW band, finite F spin couplings are 
generated which further stabilize the F-AF state, as seen from the spin wave dispersion 
in Fig. 4(a). Here A=3 and a small Hund's coupling has been included. At 
y = 8%, the spin wave energy at the F zone boundary X exceeds that at the AF 
zone boundary, as observed experimentally. The F-AF state is seen to be stable for 
very low electron doping, and gets rapidly destabilized with increasing electron doping 
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due to suppression of AF spin couplings. This behaviour is in agreement with the 
observed rapid suppression of magnetic order in iron pnictides with electron doping 
(due to F substitution of O atoms in LaOi^xFxFeAs or Ni substitution of Fe atoms in 
BaFe2_xNixAs2) [24]. 

Notably, the spin wave dispersion, energy scale, and doping behaviour for the 
electron doped SDW state are all very similar to Fig. 2 for the hole doped SDW state as 
well as the single-band model. Evidently, it is the generation of the ferromagnetic spin 
couplings in all three cases which is the common crucial factor in stabilizing the (0,7r) 
ordered SDW state. For A=2, the SDW state effective gap ~ 400meV is well above the 
maximum spin wave energy, which prevents spin wave excitations from rapidly decaying 
into a particle-hole continuum, as indeed not observed experimentally up to energies of 
200 meV. 

4. Conclusions 

Spin wave energies were obtained in the (0, tt) ordered SDW state of the minimal 
two-band model for iron pnictides including the Hund's coupling term between the 
degenerate dxz and dyz Fe orbitals. Negative spin wave energies were taken as signalling 
instability of the (0, vr) SDW state with respect to transverse spin fluctuations, indicating 
significantly weakened F or AF spin couplings. A robust (0,7r) SDW state was found 
in two distinct doping regimes — finite hole doping in the lower SDW band for small 
second neighbour hoppings (squarish electron/hole pockets), and small electron doping 
in the upper band for comparable hopping terms (circular electron/hole pockets). 

In both cases, the Hund's coupling term J was found to strongly stabilize the 
(0,7r) SDW state due to the generation of additional ferromagnetic spin couplings 
J'^[x^]'^f involving the inter-orbital part of the particle-hole propagator. Furthermore, 
the spin wave dispersion, energy scale, and doping behaviour in both cases were found 
to be very similar to that for the one-band model, which was ascribed to the carrier- 
induced ferromagnetic spin couplings as the common crucial factor in stabilizing the 
(0,7r) ordered SDW state. The emergence of F spin couplings at finite hole/electron 
doping results in a characteristic peak spin- wave energy at {11,11), in agreement with 
neutron scattering experiments on iron pnictides, and also accounts for the large planar 
anisotropy observed between the AF and F spin couplings in the ab plane. 

Evidence of Fermi surface folding associated with the SDW state has been observed 
in recent ARPES studies [23]. Electronic quasiparticle dispersion and spectral function 
renormalization in the (0, tt) ordered SDW state due to electron-magnon interaction and 
multiple magnon emission- absorption processes should therefore be of interest. 
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